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Exact  generation  of  epsilon-efficient  solutions 
in  multiple  objective  programming 


Alexander  Engau*  and  Margaret  M.  Wiecek* 
Department  of  Mathematical  Sciences 
Clemson  University,  SC,  USA 


Abstract 

It  is  a  common  characteristic  of  many  multiple  objective  programming  problems 
that  the  efficient  solution  set  can  only  be  identified  in  approximation:  since  this  set 
often  contains  an  infinite  number  of  points,  only  a  discrete  representation  can  be  com¬ 
puted,  and  due  to  numerical  difficulties,  each  of  these  points  itself  might  in  general  be 
only  approximate  to  some  efficient  point.  From  among  the  various  approximation  con¬ 
cepts,  this  paper  considers  the  notion  of  epsilon-efficiency  which  has  also  been  shown  to 
be  of  relevance  other  than  merely  for  the  purpose  to  approximate  solutions.  Following 
preceding  work  by  the  same  authors,  new  generating  methods  are  proposed  to  resolve 
various  drawbacks  of  those  methods  derived  earlier.  Supporting  theoretical  results  are 
established  and  the  methods  demonstrated  on  an  engineering  design  example. 

Keywords:  multiple  objective  programming,  epsilon-efficient  solutions,  epsilon-Pareto 
outcomes,  approximation 


1  Introduction 

Mathematical  optimization  and  decision  making  problems  with  multiple  objectives  have  been 
studied  for  more  than  fifty  years  now,  see  Koopmans  (1951),  and  many  results  have  been 
established  in  theory,  methodology  and  applications  as  collected  in  various  monographs  by 
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of  Excellence  for  Modeling  and  Simulation  of  Ground  Vehicles  at  the  University  of  Michigan,  and  by  the 
National  Science  Foundation,  grant  number  DMS-0425768. 
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Chankong  and  Haimes  (1983),  Sawaragi  et  al.  (1985),  Yu  (1985),  Steuer  (1986)  or  more 
recently  Miettinen  (1999)  and  Ehrgott  (2005).  This  paper  joins  the  effort  to  promote  the 
relevance  of  epsilon-efficient  solutions  introduced  into  multiple  objective  programming  by 
Loridan  (1984)  and  White  (1986)  and  derives  a  methodology  for  their  generation. 

The  general  goal  in  any  optimization  or  decision  making  process  is  to  identify  a  single 
or  all  best  solutions  within  a  set  of  feasible  points  or  alternatives.  While  it  is  theoretically 
possible  to  identify  the  complete  set  of  best  solutions,  finding  an  exact  description  of  this  set 
often  turns  out  to  be  practically  impossible  or  at  least  computationally  too  expensive,  and 
thus  many  research  efforts  focus  on  approximation  concepts  and  procedures,  see  Ruzika  and 
Wiecek  (2005).  Therefore,  in  practice  the  optimizing  decision  maker  is  frequently  satisfied 
with  suboptimal  solutions,  provided  the  loss  in  optimality  can  be  justified  by  significant 
gain  with  respect  to  model  simplicity  and  computational  benefits.  The  consideration  of  such 
tradeoff  reflects  the  common  belief  that  the  concept  of  epsilon-efficient  solutions  accounts 
for  modeling  limitations  or  computational  inaccuracies,  and  thus  is  tolerable  rather  than 
desirable.  Consequently,  solution  methods  purposely  avoiding  efficiency  while  seeking  to 
guarantee  epsilon-efficiency  have  not  been  well  developed. 

Following  earlier  work  by  Engau  and  Wiecek  (2005b)  who  investigate  the  significance 
of  epsilon-efficient  solutions  in  practical  decision  making  situations,  this  paper  continues 
the  development  of  methods  capable  of  generating  epsilon-efficient  solutions  for  multiple 
objective  programming  problems.  Consider  Figure  1. 


Figure  1:  Generating  methods  for  e-efficient  solutions 
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In  the  past,  much  effort  has  been  undertaken  to  develop  applicable  methods  to  solve  a 
multiple  objective  program  MOP  to  efficiency,  and  perhaps  the  most  prominent  approach 
among  those  is  to  formulate  an  auxiliary  single  objective  program  SOP  (A)  which  can  be 
solved  to  optimality  (B)  using  traditional  linear  or  nonlinear  optimization  techniques  and 
whose  optimal  solutions  give  rise  to  an  associated  efficient  solution  for  the  MOP  (C).  Based  on 
the  concepts  of  e-optimality  and  e-efficiency  for  single  and  multiple  objective  programming, 
respectively,  several  authors  also  investigated  relationships  between  e  G  M  and  e  G  Mm  for 
which  an  e-optimal  solution  for  the  SOP  (D)  produces  an  e-efficient  solution  for  the  MOP 
(E),  see  Engau  and  Wiecek  (2005b)  and  the  references  therein.  However,  these  approximate 
methods  suffer  from  the  fact  that  an  arbitrary  e-efficient  solution  for  the  SOP,  in  general, 
does  not  guarantee  a  full  relaxation  of  efficiency  for  the  MOP,  while  this  may  form  the 
intended  purpose  of  generating  such  solutions  as  mentioned  above.  In  order  to  resolve  this 
drawback,  we  propose  alternative  methods  for  the  exact  generation  of  those  fully  relaxed  e- 
efficient  solutions  for  the  MOP,  and  more  precisely,  we  introduce  four  possible  modifications 
of  the  original  MOP  into  an  £-MOP  (F)  so  that  optimal  solutions  (H)  for  its  associated 
scalarized  problem  (G)  yield  associated  ^-efficient  solutions  (I)  to  the  original  MOP.  Based 
on  a  thorough  literature  review,  we  believe  that  this  approach  and  its  relevance  have  not 
been  studied  before  and  thus  are  new  to  the  research  community. 

The  organization  of  the  remaining  text  is  now  as  follows.  The  adopted  notation  and  basic 
terminology  with  special  emphasis  on  the  concepts  of  efficiency  and  ^-efficiency  for  the  MOP 
and  e-optimality  for  its  scalarized  counterpart  are  established  in  Section  2  that  also  revisits 
and  summarizes  some  prior  results  from  Engau  and  Wiecek  (2005b).  Section  3  deals  with 
the  development  of  exact  generating  methods  and  gives  further  motivation  and  supporting 
theoretical  results.  In  Section  4,  these  new  methods  are  applied  to  an  engineering  design 
problem  and  compared  to  the  previous  approach  in  order  to  demonstrate  the  significant 
improvements  achieved.  Final  remarks  and  possible  future  research  directions  are  provided 
in  the  concluding  Section  5. 

2  Notation,  basic  terminology  and  prior  results 

We  let  fi  :  Mn  — >  R,  i  =  1  be  m  real-valued  objective  functions  and  define  the 

vector-valued  function  /  :  Mn  — >  by  f(x)  :=  [fi(x), . . . ,  fm(x)}T .  Then  we  formulate  the 
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minimization  problem 


MOP:  minimize  f(x)  subject  to  x  G  X  C  Mn 

as  the  multiple  objective  program  of  interest,  for  which  X  is  a  given  set  of  feasible  alternatives 
in  the  decision  space  Mn.  A  feasible  decision  x  G  X  is  evaluated  by  the  objective  function  / 
to  produce  the  outcome  f{x)  in  the  objective  or  outcome  space  Mm.  We  define  the  set  of  all 
attainable  outcomes  as  the  image  of  X  under  /  and  denote  this  set  by 

Y  :=  /( X)  :=  {y  G  :  y  =  f(x)  for  some  x  G  X}. 

2.1  Efficient  and  ^-efficient  solutions 

We  first  define  the  underlying  concept  of  minimality  on  Mm.  Given  two  outcomes  y 1  and 
y2  G  Y  C  Mm,  we  write  y1  A  y2  if  and  only  if  yj  <  y2  for  all  i  —  1, . . . ,  m,  y1  <  y2  if  and 
only  if  yl  A  y2  and  y1  ^  y2,  y1  <  y2  if  and  only  if  yj  <  y2  for  all  i  —  1, . . . ,  m,  and  y1  A  y2, 
y1  >  y2,  y1  >  y2  if  and  only  if  y2  ^  y1,  y2  <  y1,  y2  <  y1,  respectively.  Then  x°  G  X  is 
called  an  efficient  or  weakly  efficient  decision  for  the  MOP  if  and  only  if  there  does  not  exist 
another  feasible  x  G  X  such  that  f(x)  <  f(x° )  or  f(x)  <  f(x°),  and  the  sets  of  all  efficient 
or  weakly  efficient  solutions  are  denoted  by  X°  or  X?° ,  respectively.  An  efficient  or  weakly 
efficient  solution  x°  G  X°  or  produces  a  Pareto  or  weak  Pareto  outcome  y°  =  f(x°),  and 
the  sets  of  all  Pareto  or  weakly  Pareto  outcomes  are  denoted  by  Y°  or  Y°,  respectively. 

Solving  the  MOP  now  means  to  find  the  set  (or  a  subset)  of  (weakly)  efficient  decisions 
or  (weak)  Pareto  outcomes.  Conditions  for  the  existence  of  these  solutions  are  established 
by  Hartley  (1978),  Corley  (1980),  Borwein  (1983)  and  Sawaragi  et  al.  (1985)  and,  in  general, 
can  be  guaranteed  under  certain  compactness  assumptions  on  the  outcome  set  Y.  For  a 
comparison  of  several  existence  results,  we  refer  to  the  recent  survey  provided  by  Sonntag 
and  Zalinescu  (2000). 

Although  theoretical  studies  have  derived  convenient  representations  for  the  set  of  Pareto 
outcomes  in  the  objective  space,  it  is  important  to  realize  that  given  a  Pareto  outcome 
y°  G  Y°  or  any  other  attainable  outcome  y  G  Y,  it  is  in  general  not  possible  to  easily 
identify  an  efficient  decision  x°  G  X°  or  feasible  decision  x  G  X  that  produces  this  outcome 
y°  —  f(x°)  or  y  —  f(x),  respectively. 

Given  £  G  Mm,  £  2  0,  a  feasible  decision  x°  G  X  is  called  an  e-efficient  or  weakly  e- 
efficient  solution  for  the  MOP  if  and  only  if  there  does  not  exist  another  feasible  x  G  X  such 
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that  f(x)  <  f(x° )  —  £  or  f(x)  <  f(x°)  —  e,  and  the  sets  of  all  e-efficient  or  weakly  e-efficient 
solutions  are  denoted  by  or  Xw  \  respectively.  An  e-efficient  or  weakly  e-efficient  solution 
x°  G  X&  or  xif  '1  produces  an  e-Pareto  or  weak  e-Pareto  outcome  y°  =  f(x°),  and  the  sets 
of  all  e-Pareto  or  weakly  e-Pareto  outcomes  are  denoted  by  or  Yw  \  respectively.  In 
particular,  we  obtain  that  =  X°,  Xw  ^  =  X =  Y°  and  =  Y°. 


Figure  2:  Sets  of  (weakly)  Pareto  and  (weakly)  e-Pareto  outcomes 


Figure  2  illustrates  the  relationships  between  the  different  solutions  concepts  for  a  generic 
outcome  set  of  some  biobjective  problem.  Consider  the  depicted  set  7  C  l2.  The  set  Y°  of 
Pareto  outcomes  is  given  by  the  curve  connecting  points  A  and  B,  while  all  weakly  Pareto 
outcomes  of  the  set  Y£  lie  on  the  extended  curve  from  C  to  D.  Adding  the  given  vector 
£  G  M2  now  shifts  the  points  A,  B,  C  and  D  to  E,  F,  G  and  H,  respectively,  and  then  the 
set  y(£)  of  £-Pareto  outcomes  is  the  shaded  area  enclosed  by  all  marked  points  without  the 
curve  connecting  I,  G,  E,  F,  H  and  J,  while  the  set  Yyf'1  of  weak  £-Pareto  outcomes  includes 
the  curve  from  1  to  J.  It  follows  that  for  this  particular  example,  the  set  Y ^  is  closed, 
while  the  set  of  £-Pareto  outcomes  is  neither  closed  nor  open.  Consequences  of  this  last 
observation  are  further  addressed  at  the  beginning  of  Section  3. 


2.2  Scalarizing  the  MOP  and  e-optimal  solutions 

The  most  frequent  approach  to  solving  MOPs  is  to  use  a  parametric  characterization  of  the 
set  of  efficient  decisions  by  aggregating  the  objective  functions  fi,  i  =  1  into  one 
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real-valued  objective  s  which  is  minimized  if  and  only  if  the  associated  decision  x  G  X  is 
(weakly)  efficient  for  the  original  MOP.  We  call  such  an  objective  function  s  a  scalarization 
function  for  the  MOP  and  require  that  this  function  preserves  the  partial  order  defined  on 
Y  C  Mm. 

Let  y 1  and  y2  G  Y  be  any  two  attainable  outcomes.  Then  we  call  the  scalarization 
function  s  increasing  if  y 1  ^  y2  and  y1  <  y2  implies  that  s(y1)  <  s(y2)  and  s^1)  <  s(y2), 
respectively.  If.  in  addition,  y1  <  y2  implies  that  s^y1)  <  s(y2),  then  s  is  called  strictly 
increasing. 

Based  on  several  ideas  and  different  methodological  concepts,  a  large  number  of  increasing 
scalarization  functions  is  available,  possibly  depending  on  additional  scalarization  parame¬ 
ters  or  auxiliary  variables.  Common  characteristics  of  these  functions  include  the  use  of 
weighting  coefficients  as  in  the  weighted  sum  approach  by  Geoffrion  (1968),  reference  points 
as  in  Benson  (1978)  or  in  combination  with  norms  such  as  the  weighted  and  augmented  ip 
norms  as  in  Stener  (1986),  objective  function  levels  as  in  the  constrained  objective  approach 
discussed  in  Chankong  and  Haimes  (1983),  search  directions  as  in  Pascoletti  and  Serafim 
(1984),  aspiration  levels  as  in  goal  programming  approaches,  or  other  combinations  of  the 
above.  Wierzbicki  (1986)  examines  several  properties  of  scalarization  functions,  including 
their  ability  to  generate  the  complete  set  of  efficient  decisions  when  varying  the  respective 
scalarization  parameters.  Li  et  al.  (1999)  derive  relationships  between  these  parameters  for 
four  different  scalarization  methods  when  generating  the  same  efficient  solution. 

For  the  purpose  of  this  paper,  we  formally  define  a  scalarization  function  as  an  increasing 
function  s  :  Y  x  II  — >  R,  where  II  is  the  set  of  admissible  scalarization  parameters.  Then  the 
parametric  single  objective  program  associated  with  the  MOP  is  the  minimization  problem 

minimize  s(y,i r)  subject  to  y  G  Y 

with  7r  G  II  as  fixed  scalarization  parameter.  However,  since  the  set  Y  is  in  general  not  given 
in  explicit  form,  the  minimization  needs  to  be  done  over  the  feasible  set  X  for  the  MOP, 
thus  yielding  the  final  formulation 

SOP(7r):  minimize  s(f(x),  n)  subject  to  x  G  X. 

Then  x°  G  X  and  y°  =  f(x°)  G  Y  are  called  optimal  solutions  for  the  SOP  if  and  only  if 
s(f(x°),  7r)  <  s(f(x ),  7r )  for  all  x  G  X.  Given  e  G  R,  e  >  0,  then  x°  G  X  and  y°  =  f(x° )  G  Y 
are  called  e-optimal  or  strictly  e-optimal  solutions  for  the  SOP  if  and  only  if  s(f(x°),  n)  < 
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s(f(x),  it)  +  e  or  s(f(x°),  it)  <  s(f(x),  7 r)  +  e  for  all  x  G  X,  respectively. 

Note  that  we  define  e-efficient  and  weakly  e-efficient  solutions  for  the  MOP,  but  e-optimal 
and  strictly  e-optimal  solutions  for  the  SOP,  in  accordance  with  the  standard  terminology 
established  in  the  literature,  here  namely  Loridan  (1982,  1984)  and  White  (1986). 

The  following  lemma  illustrates  the  importance  for  the  chosen  scalarization  function  to 
be  increasing  in  order  to  relate  optimal  solutions  for  the  SOP  and  efficient  decisions  for  the 
MOP  and  is  similarly  established  by  Wierzbicki  (1986). 

Lemma  2.1.  Let  the  MOP  and  an  SOP  with  increasing  scalarization  function  s  be  given. 
If  x°  €  X  is  optimal  for  the  SOP,  then  x°  is  weakly  efficient  for  the  MOP.  If  s  is  strictly 
increasing,  then  x°  is  efficient. 

For  different  choices  of  the  scalarization  parameter  7r  €  n,  in  general  we  obtain  different 
optimal  solutions  for  the  SOP  and  thus  different  efficient  decisions  for  the  MOP. 

Two  of  the  most  popular  scalarization  methods  are  the  weighted-sum  and  the  weighted- 
Tchebycheff  norm  scalarization  which  we  formulate  now  for  later  reference  in  the  subsequent 
section. 

Given  the  MOP,  let  w  G  Mm,  w  >  0  be  a  given  vector  of  weights  for  the  objective  function 
/  =  (fi, . . . ,  fm)T.  Then  the  weighted-sum  scalarization  for  the  MOP  is  defined  as 

m 

WS(w):  minimize  Wjfj(x)  subject  to  x  G  X. 

1=1 

In  addition  to  the  above,  let  r  G  be  a  given  reference  point.  Then  the  weighted- 
Tchebycheff  norm  scalarization  for  the  MOP  is  defined  as 

TN(r,  w):  minimize  max  {u!t(ft(x)  —  r,)}  subject  to  x  G  X. 

i=l,...,m 

2.3  Approximate  generating  methods 

As  we  mentioned  before,  it  is  theoretically  possible  to  identify  the  complete  set  of  efficient 
decisions.  In  practice,  however,  often  only  a  discrete  representation  of  the  efficient  set  can 
actually  be  computed  and,  due  to  numerical  difficulties,  each  of  these  points  itself  might 
only  be  approximate  to  some  truly  efficient  point.  In  Engau  and  Wiecek  (2005b),  seven 
different  scalarization  methods  are  examined  with  respect  to  e-optimality  and  e-efficiency  and 
corresponding  relationships  are  established  between  e  G  M  and  e  G  Mm  when  solving  both 
the  scalarized  problem  SOP  for  an  e-optimal  and  the  original  MOP  for  an  e-efficient  solution. 
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Based  on  these  results  a  generic  procedure  for  the  computation  of  ^-efficient  decisions  for 
the  MOP  is  proposed,  thereby  relying  on  the  idea  that  solving  any  of  these  SOPs  with  some 
scalarization  parameter  n  G  fl  for  a  (strictly)  e-optimal  solution  also  produces  a  (weakly) 
e-efficient  decision  for  the  MOP,  according  to  the  path  ADE  in  Figure  1. 

This  procedure  assumes  that  based  on  prior  knowledge  and  experience,  a  decision  maker 
is  able  to  choose  both  the  vector  parameter  e  and  a  suitable  SOP  together  with  the  required 
scalarization  parameter  n.  Clearly,  these  choices  depend  on  the  actual  type  of  problem  to 
be  solved  and  the  desired  relaxation  of  efficiency  for  the  MOP.  Furthermore,  all  results  only 
imply  an  upper  bound  on  the  set  of  admissible  values  for  e  such  that  an  e-optimal  solution  for 
the  SOP  is  guaranteed  to  be  (weakly)  e-efficient  for  the  MOP,  and  hence,  when  solving  the 
SOP  for  an  e-optimal  solution,  it  is  recommended  that  the  actual  e  should  be  chosen  equal 
to  the  upper  bound  if  the  full  relaxation  is  to  be  accomplished.  However,  while  the  decision 
maker  may  gain  certain  control  over  the  e-efficiency  of  generated  feasible  decisions  for  the 
MOP  by  varying  the  parameter  e,  this  approach  highly  depends  on  the  capability  to  find 
reasonably  relaxed  e-optimal  solutions  for  the  SOP.  Although  in  principle  accomplishable  by 
a  simple  modification  of  the  termination  criterion  of  the  optimization  algorithms  used,  the 
lack  of  accurate  bounds  on  the  actual  optimum  and  thus  the  uncertainty  regarding  the  actual 
relaxation  achieved  motivates  the  development  of  methods  for  which  e-efficient  decisions  can 
be  obtained  as  (exact)  efficient  decisions,  thereby  guaranteeing  the  full  desired  relaxation. 

3  Exact  generating  methods 

We  formulate  the  particular  problem  of  interest  as  to  find  the  worst  solutions  among  the  set  of 
e-efficient  decisions  and  then  introduce  four  different  approaches  to  modify  the  original  MOP 
formulation  so  that  an  exact  solution  to  the  modified  problem  is  such  fully  relaxed  e-efficient 
decision  for  the  original  MOP.  For  obvious  reasons  and  in  order  to  avoid  pathological  results, 
we  exclude  the  case  e  =  0  from  all  the  following  discussion.  Furthermore,  we  mentioned 
before  that  the  existence  of  efficient  solutions  can,  in  general,  be  guaranteed  for  a  closed  and 
bounded  set  of  feasible  points  in  the  outcome  space,  while  the  discussion  of  Figure  2  revealed 
that  the  set  of  e-Pareto  outcomes  does  not  necessarily  satisfy  these  conditions.  Hence,  in 
order  to  avoid  emptiness  of  the  set  of  efficient  solutions  for  the  problem  to  be  formulated, 
we  choose  the  set  of  weakly  e-efficicnt  decisions  as  the  underlying  feasible  set. 


Let  the  MOP  and  e  G  Mm,  e  >  0  be  given.  The  associated  e-MOP  is  defined  as 
e-MOP:  maximize  f(x)  subject  to  x  G  X^\ 

Analogously  to  the  definition  of  efficient  decisions  for  the  minimization  MOP,  a  feasible 
decision  x°  G  X if'1  is  called  an  efficient  or  weakly  efficient  solution  for  the  e-MOP  if  and 
only  if  there  does  not  exist  another  feasible  x  G  Xw ^  such  that  f(x°)  <  f{x)  or  f(x°)  < 
f(x),  respectively.  The  sets  of  all  efficient  or  weakly  efficient  solutions  for  the  e-MOP  are 
denoted  by  e-X°  and  e-X °,  respectively,  and  analogously  we  define  e-Y°  and  e-Y°  as  the 
corresponding  Pareto  and  weak  Pareto  sets. 

Note  that  the  set  of  weakly  e-efficient  decisions  serving  as  the  feasible  decision  set  for 
the  e-MOP  is  in  general  not  known,  and  clearly,  without  knowledge  of  this  set,  solving  the 
e-MOP  is  in  general  not  possible.  In  another  paper,  Engau  and  Wiecek  (2005a)  study  the 
problem  of  maximizing  over  the  set  of  e-efficient  solutions  in  a  more  abstract  setting  and 
propose  an  alternative  characterization  of  the  efficient  solution  set  for  the  e-MOP.  To  clarify 
terminology,  we  first  define  e-solutions  as  fully  relaxed  e-Pareto  outcomes  in  the  following 
sense. 

Let  the  MOP  and  e  G  Mm,  e  >  0  be  given.  Then  an  outcome  y*  G  Y  is  called  an 
e-outcome  if  there  exists  a  weak  Pareto  outcome  y°  G  Y°  such  that  y*  =  y°  +  e,  and  the  set 
of  all  e-outcomes  is  denoted  by  Y*.  A  decision  x*  G  X  is  called  an  e-decision  if  it  produces 
an  e-outcome,  y*  =  f(x*)  G  Y*,  and  the  set  of  all  e-decisions  is  denoted  by  A"*.  Formally, 
we  write 

Y*  =  Y  O  (Y°  +  e)  and  X*  =  {iG  X  :  f(x)  G  Y*} 

We  illustrate  the  new  solutions  concepts  and  sets  for  the  generic  outcome  set  given  in 
Figure  2.  Recall  that  the  set  of  Pareto  outcomes  Y°  is  given  by  the  curve  connecting  points 
A  and  B,  while  all  weakly  Pareto  outcomes  in  the  set  Y°  lie  on  the  curve  from  C  to  D.  Adding 
the  given  vector  £Gt2  then  yields  the  set  Y*  of  e-outcomes  as  the  curve  connecting  points 
G  and  H.  As  before,  the  set  Yw  ^  of  weakly  e-Pareto  outcomes  is  the  shaded  area  enclosed  by 
all  marked  points  and  now  forms  the  closed  set  of  attainable  outcomes  for  the  e-MOP.  Then 
the  Pareto  set  for  the  e-MOP  is  given  by  the  curve  from  E  to  F  without  the  two  endpoints 
E  and  F,  but  including  the  two  points  I  and  J,  and  finally  the  set  of  weak  Pareto  outcomes 
e-Y°  consists  of  all  points  along  the  complete  curve  connecting  I  and  J. 

Moreover,  Figure  2  now  suggests  that  the  set  of  weakly  efficient  solutions  e-X°  for  the  e- 
MOP  contains  the  set  of  those  weakly  e-efficient  decisions  for  the  MOP  that  are  not  e-efficient 
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and,  moreover,  that  this  set  includes  the  set  of  e-decisions  for  the  MOP  as  a  (possibly  proper) 
subset.  This  intuition  is  confirmed  as  special  case  of  the  results  in  the  aforementioned  paper. 

Lemma  3.1.  Let  the  MOP  and  £  G  Mm ,  e  >  0  be  given.  Then  all  e-decisions  and  e-outcomes 
for  the  MOP  belong  to  the  sets  of  weakly  efficient  decisions  and  weak  Pareto  outcomes  for 
the  associated  e-MOP, 

X *  C  X$  \  X(£)  C  £-X°w  and  Y*  C  Y,^  \  F(e)  C  e-Y°. 

Motivated  by  this  lemma,  for  the  remaining  discussion  we  continue  to  focus  on  weak 
efficiency  as  the  underlying  solution  concept  and,  in  order  to  solve  the  e-MOP,  now  pursue 
the  development  of  methods  for  finding  e-decisions  and  e-outcomes  for  the  original  MOP. 
A  steady  assumption  in  each  of  the  following  four  sections  is  that  it  is  always  possible  to 
generate  one  or  more  initial  (weakly)  efficient  decisions  x°  G  X°,  producing  a  (weak)  Pareto 
outcome  y°  =  f(x°)  G  Y°  for  which  y*  =  y°  +  e  G  Y*  is  an  e-outcome  in  the  set  Y  of 
attainable  outcomes.  Then,  by  Lemma  3.1,  this  outcome  y*  G  Y*  is  also  (weak)  Pareto  for 
the  e-MOP,  y*  G  e-Y°. 

We  emphasize  that  the  main  difficulty  is  not  to  fold  an  e-outcome  y*  in  the  set  of  at¬ 
tainable  outcomes,  but  to  identify  an  associated  e-decision  x*  G  A"*  that  produces  this 
e-outcome  y*  =  f(x*)  G  Y*  C  e-Y°.  Now  the  following  four  approaches  are  proposed  as 
possible  solution  methods. 

3.1  Constrained  objective  approach 

Recall  the  formulation  of  the  MOP  and  the  e-MOP 

MOP:  minimize  f(x)  subject  to  x  G  X  C  Mn; 
e-MOP:  maximize  f(x)  subject  to  x  G  X^\ 

Now  let  y°  G  Y°  be  a  (weak)  Pareto  outcome  for  which  the  associated  e-outcome  y*  = 
y°  +  £  E  Y*  is  attainable  and  consider  the  multiple  objective  problem 

CO (y*):  maximize  f(x)  subject  to  f(x)  ^  y*,x  G  X. 

Proposition  3.2.  All  efficient  solutions  for  the  CO  (y*)  are  e- decisions  for  the  original  MOP. 
Proof.  Observe  that  CO (y*)  has  the  point  y*  as  a  unique  Pareto  outcome,  which  by  as- 
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sumption  is  an  attainable  ^-outcome,  y*  =  y°  +  e  G  Y*.  Thus  all  efficient  decisions  x*  for 
the  CO (y*)  must  satisfy  that  f(x*)  =  y*,  and  hence  all  efficient  decisions  x*  for  the  CO (y*) 
are  e-decisions  for  the  original  MOP,  x*  G  X* .  □ 

Hence,  Proposition  3.2  now  allows  to  solve  the  CO (y*)  in  order  to  identify  corresponding 
e-decisions  x*  G  A"*  for  all  e-outcomes  y*  G  Y*,  which,  by  Lemma  3.1,  belong  to  the  (weakly) 
efficient  solution  set  for  the  e-MOP,  x*  G  e-X°,  and  y*  G  s-Y°. 

Note  that  the  above  problem  formulation  resembles  the  constrained  objective  scalariza- 
tion 

minimize  f\  (x)  subject  to  fi(x)  <  y*,  i  —  2, . . . ,  m,  x  G  X 

for  which  all  but  one  objectives  are  bounded  from  above  to  guarantee  that  the  objective 
function  values  improve  or  at  least  satisfy  certain  aspiration  levels.  However,  in  the  formu¬ 
lation  of  the  CO (y*)  all  original  objectives  remain  in  the  multiple  objective  function  while 
the  objective  function  levels  implied  by  y*  take  the  role  of  reservation  levels  which  must  not 
be  exceeded.  Finally,  observe  that  once  the  problem  is  formulated,  it  can  be  solved  using 
any  appropriate  scalarization  method. 

3.2  Reference  point  approach 

Perhaps  the  main  difficulty  for  finding  ^-decisions  is  that  the  associated  e-outcomes,  in 
general,  are  interior  points  of  the  outcome  set,  so  that  many  optimization  procedures  that 
depend  on  the  fact  that  optimal  solutions  occur  at  the  boundary  cannot  be  used.  Therefore 
the  above  constrained  objective  approach  restricts  the  set  of  attainable  outcomes  in  such  a 
way  that  an  e-outcome  can  be  found  at  the  boundary  of  the  modified  outcome  set.  As  a 
second  possibility  we  suggest  to  make  use  of  the  reference  point  methodology  in  multiple 
objective  programming  for  which,  very  recently,  Lin  (2005)  investigates  the  consequences 
of  the  location  of  the  chosen  reference  point  with  respect  to  the  efficiency  of  the  solutions 
generated.  We  base  our  approach  on  the  weighted  £p  norm 

\  !/P 

Wi\yi  -  n\p 

which  uses  a  reference  point  r  G  M"1  and  a  weighting  parameter  w  G  Mm,  w  >  0.  Here 
r  G  is  usually  chosen  to  be  the  ideal  or  a  utopia  point,  and  then  the  absolute  values  are 
usually  dropped  since  r  A  y  for  all  y  G  Y.  In  particular  for  p  =  oo,  the  weighted  norm 
is  defined  to  be  the  weighted-Tchebycheff  norm  introduced  in  Section  2.2. 


tP{y,r,w)  = 


m 

E 

,  i= 1 


li 


Now  let  y°  G  Y°  be  a  (weak)  Pareto  outcome  for  which  the  associated  e-outcome  y*  = 
y°  +  e  G  Y*  is  attainable  and  choose  positive  weights  w  >  0  and  r  =  y*  as  reference  point 
for  the  single  objective  reference  point  problem 


RP  (y*  ,w):  minimize 


i/p 


Wi\fi 


X 


Vi 


*\p 


subject  to  x  G  X. 


Proposition  3.3.  All  optimal  solutions  for  the  R P(y*,w)  are  e-decisions  for  the  original 

MOP. 


Proof.  Observe  that  since  w  >  0,  RP (y*,w)  has  the  point  y*  as  a  unique  optimal  solution, 
which  by  assumption  is  an  attainable  e-outcome,  y*  =  y°+e  G  Y*.  Thus  all  optimal  solutions 
x*  for  the  RP(i/*,tc)  must  satisfy  that  f(x*)  =  y*,  and  hence  all  optimal  solutions  x*  for  the 
RP(2/*,w)  are  e-decisions  for  the  original  MOP,  x*  G  A"*.  □ 

Hence,  Proposition  3.3  now  allows  to  solve  the  RP(y*,w)  in  order  to  identify  correspond¬ 
ing  e-decisions  x*  G  A"*  for  all  e-outcomes  y*  G  Y*,  which,  by  Lemma  3.1,  belong  to  the 
(weakly)  efficient  solution  set  for  the  e-MOP,  x*  G  e-X“  and  y*  G  e-Y°. 


3.3  Projection  approach 

As  mentioned  earlier,  most  traditional  optimization  methods  fail  to  identify  points  in  the 
interior  of  the  outcome  set  and  thus  need  to  be  modified  to  allow  for  the  generation  of  e- 
decisions  or  e-outcomes  which,  in  general,  do  not  occur  at  the  boundary.  In  this  section 
we  investigate  a  preliminary  projection  approach  which,  besides  offering  theoretical  insight, 
prepares  for  more  elaborate  methods  in  the  subsequent  and  final  section.  The  main  idea  is 
to  introduce  an  additional  objective  fm+ 1  :  IP”  — >  IP  such  that  interior  points  of  the  original 
outcome  set  Y  C  become  boundary  points  of  an  augmented  outcome  set  y  C  Mm+1. 
The  idea  is  illustrated  by  the  following  example. 

Example  3.4.  Let  the  MOP  with  objective  function  /  =  (/1; . . . ,  /m)r  be  given,  and  intro¬ 
duce  an  additional  component  by  setting  fm+1(x)  =  0  for  all  x  G  A".  Then 

y  =  Y  x  {0}  =  f(X)  x  {0}  C  Mm+1 

is  contained  in  a  hyperplane  in  m  +  1  dimensional  space,  and  hence  all  outcomes  y  G  Y 
become  boundary  points  (y,  0)r  G  y  in  the  augmented  outcome  space  IPm+1. 
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Now  let  y°  G  Y°  be  a  (weak)  Pareto  outcome  for  which  the  associated  e-outcome  y*  = 
y°  +  £  G  Y*  is  attainable  and  consider  the  single  objective  problem 


PP (y*,p):  minimize  p  subject  to 


(p)  +  '*  (-1) €  y  =  f(x)  x  {0} 


with  arbitrary  p  G  M.  Observe  that  solving  this  problem  is  equivalent  to  finding  the  projec¬ 
tion  of  the  point  (y*,  p)T  G  x  ffi.  onto  the  set  y  and  has  p  =  p  as  solution  with  associated 
outcome  ( y *,  0)  G  y. 

Proposition  3.5.  All  optimal  solutions  for  the  PP (y*,p)  yield  e-decisions  for  the  original 

MOP. 


Proof.  Observe  that  p  —  p  is  the  unique  feasible  and  thus  optimal  value  for  the  PP (y*,p). 
The  associated  solution  is  ( y *,  0)T  G  f(X )  x  {0}  and  thus  will  be  produced  by  an  associated 
e-decision  x*  for  which  f(x*)  =  y* .  □ 

While  clearly  trivial  from  an  analytical  point  of  view,  Proposition  3.5  now  allows  to  solve 
the  PP(y*,p)  numerically  in  order  to  identify  corresponding  e-decisions  x  G  X*  for  all  e- 
outcomes  y*  G  Y*,  which,  by  Lemma  3.1,  belong  to  the  (weakly)  efficient  solution  set  for  the 
e-MOP,  x*  G  s-X°  and  y*  G  e-Y°. 

The  practical  applicability  of  the  above  problem  is  not  immediate,  however,  since  the 
optimal  solution  is  in  particular  the  only  feasibly  solution.  Nevertheless,  the  idea  of  aug¬ 
menting  the  outcome  space  and  simultaneously  lifting  the  set  of  attainable  outcomes  seems 
quite  appealing  and  thus  motivates  the  following  discussion  of  a  general  set  lifting  approach 
as  one  further  method  for  the  generation  of  e-efficient  solutions. 


3.4  Set  lifting  approach 

We  continue  to  investigate  the  augmentation  idea  described  in  the  previous  section  and  pro¬ 
pose  the  formulation  of  an  augmented  MOP  (AMOP).  Then  we  study  the  AMOP  with  re¬ 
spect  to  two  scalarization  methods,  namely  the  weighted-sum  and  the  weighted-Tchebycheff 
norm  scalarization,  and  show  how  ^-decisions  for  the  MOP  can  be  identified  as  efficient 
decisions  for  the  respective  AMOPs. 

Let  the  MOP  be  given,  and  let  a  :  Y  C  Mm  ->Rbea  decreasing  scalarization  function, 
a(yl )  >  a(y2)  whenever  yl,y2  G  Y,  y 1  ^  y2 . 
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Then  the  augmented  multiobjective  program  is  defined  as 


AMOP:  minimize  [f(x),  a(f(x))}T  subject  to  x  G  X, 

and  the  function  a  is  called  the  augmentation  function  of  the  AMOP. 

The  next  lemma  shows  how  increasing  scalarization  functions  can  be  used  to  define  suit¬ 
able  augmentation  functions  for  the  AMOP  and  prepares  the  two  subsequent  main  theorems 
of  this  section. 

Lemma  3.6.  Let  s  :  Y  C  — >  R  be  an  increasing  scalarization  function,  and  let  y°  G  Y 

be  an  optimal  solution  for  the  associated  SOP.  Define 

a(y)  ■=  (s(y)  -  sly0))-1 

and  allow  the  function  value  a(y°)  =  oo.  Then  the  function  a  :  Y  C  Mm  ^RU  {00}  defines 
a  decreasing  augmentation  function. 

Proof.  Let  yl,y2  G  Y  with  y1  A  y2  be  given.  Since  s  is  an  increasing  scalarization  function, 
this  implies  sly1)  <  s(y2)  and  thus  0  <  sly1)  —  s(y°)  <  sly2)  —  s(y°)  as  y°  is  an  optimal 
solution  for  the  SOP.  Taking  inverses  gives  (s(y x)  —  sly°))-1  >  (s(y2)  —  s(|/0))_1  and  thus 
implies  afy1)  >  afy2) .  □ 

The  two  concluding  theorems  show  that  by  solving  the  AMOP  with  the  augmentation 
function  defined  from  either  the  weighted-sum  or  the  weighted-Tchebycheff  norm  scalariza¬ 
tion,  e-decisions  for  the  MOP  can  be  obtained  as  optimal  solutions  to  the  scalarized  AMOP. 
Compare  the  illustration  of  the  weighted-sum  set  lifting  in  Figure  3  for  further  motivation 
of  this  approach. 

Theorem  3.7.  Let  the  MOP  and  e  G  e  >  0  be  given,  and  let  x°  G  X°  be  a  (weakly) 
efficient  decision  and  y°  =  f(x°)  G  Y°  be  a  (weak)  Pareto  outcome  for  the  MOP  obtained 
as  optimal  solutions  for  the  weighted- sum  scalarization  WS(tc)  with  weighting  parameter 
w  G  w  >  0.  Denote 

m 

slflx),w)  :=  ^wjilx) 

1=1 

and  let  a(y)  :=  (s(y,  w)  —  sly°,  w))_1  be  the  augmentation  function  for  the  AMOP.  Let 

m 

e  :=  s(e,w)  =  ^  WjEj. 

i=  1 
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Figure  3:  Weighted-sum  set  lifting  method 


If  the  e-outcome  y*  =  y°  +  £  G  Y*  is  attainable  for  some  feasible  x*  G  X,  y*  =  f(x*),  then 
the  e-decision  x*  G  X*  for  the  original  MOP  is  optimal  for  the  weighted-sum  scalarization 
WS^Jw,  e2]Tj  of  the  associated  AMOP. 

Proof.  Let  y*  =  y°  +  e  G  Y*  be  attainable  for  the  original  MOP  and  let  x*  G  A"*  with 
f(x*)  =  y*  be  an  associated  e-decision.  Now  consider  the  AMOP  and  the  related  weighted- 
sum  scalarization 


AMOP:  min  [f(x),(s(f(x),w)-s(f(x°),w))  ^ 

x£X 

WS([w,e2]T):  min  s(f(x),  w)  +  e2  (s(f(x),  w)  -  s(f(x°),  w))'1 , 

and  substitute  t  =  s(f(x),w)  and  t°  =  s(f(x°),w).  Then  the  objective  function  of  the 
WS  ( [w,  e2]T^\  becomes  t+e2  (t  —  f°)_1,  which  takes  its  minimum  for  t*  =  t°+e  =  s(f(x°),w)+ 
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e,  or 


m  m 

t*  =  s(y°,  w)  +  e  =  ^  wiV°i  +  Wi£i 

i— 1  i— 1 


=  Wi(yi  +  £i^>  =  s(y°  +  £>w)  =  s(y*,w)- 

i=l 

Hence  we  obtain  that  y*  G  Y*  and  thus  x*  G  A"*  are  optimal  solutions  for  the  WS^[w,  e2}T 

□ 


Other  than  the  previous  result  using  the  weighted-sum  scalarization,  the  final  theorem  is 
based  on  the  weighted-Tchebycheff  norm  and  formulated  for  any  (weakly)  efficient  decision 
x°  G  X°  or  (weak)  Pareto  outcome  y°  eY°. 

Theorem  3.8.  Let  the  MOP  and  £  G  M.m,£  >  0  be  given,  and  let  x°  G  X°  be  any  (weakly) 
efficient  decision  and  y°  =  f(x° )  G  Yf  be  the  associated  (weak)  Pareto  outcome.  Denote 

s(y,r,w):=  max  {wi(yi  -  r*)} 

mid  let  a(y)  :=  ( s(y,r,w )  —  s(y°,  r,  w))^1  be  the  augmentation  function  for  the  associated 
AMOP.  Choose  r  =  y°  as  the  reference  point  and  set 

Wi  =  ef1  for  i  —  1, . . . ,  m 

with  Wj.  =  oo  if  £i  =  0.  If  the  e-outcome  y*  =  y°  +  £  G  Y*  is  attainable  for  some  feasible 
e-decision  x*  G  X*  for  the  original  MOP  withy*  =  f(x*),  thenx*  is  optimal  for  the  weighted- 
Tchebycheff  norm  scalarization  TN([t/°,  0]t,  [w,  1]t)  of  the  AMOP. 

Proof.  Let  y*  —  y°  +  £  G  Y*  be  attainable  for  the  MOP  and  x*  G  X*  with  f(x*)  =  y*  be  a 
corresponding  e-decision.  Since  the  reference  point  for  the  augmentation  function  is  chosen 
as  r  =  y°,  we  first  obtain  that 

a(y)  =  y°,w)  -  s(yfyf  w ))_1  =  s(y,  y°,  w )-1, 
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and  hence  the  AMOP  and  the  related  weighted-Tchebycheff  norm  scalarization  are  given  by 

AMOP:  min  \f(x),  s  (f(x),  y°, 

xEX 

TN([j/°,  0]t,  [w,  1]t)  :  min  max  {wi(fi(x)  -  y°),  s  (f(x),  y°,  w)-1} 

x€X  i=l 

Now  observe  that  the  later  objective  function  can  equivalently  be  written  as 

max  max  {wi(fi(x)  -  y°)},  s  (f(x),y° ,w)_1  \ 

I  i=l,...,m  I 

and  substitute  t  =  s(f(x),y°,w )  =  maxi=i ,...,m{wi(fi(x)  —  y°)}  >  0.  Then  the  objective 
function  of  the  TN([j/°,  0]T,  [  w,  1]T)  becomes  max{t,  t  with  t  >  0,  which  takes  its  minimum 
for  t  =  1,  or 

t  =  s(y,  y°,w)=  max  {w^y,  -  y°)}  =  1. 

t=l,..., m 

Note  that  in  particular  for  y*  =  y°  +  £,  we  obtain 

s(y*,y°,  w)  =  max  {wi(y*  -  y°)}  =  max  {w^}  =  1 

t=l,...  ,m  i=l,..., m 

by  our  choice  of  wy  =  s^1,  and  hence,  y*  =  y°  +  e  e  Y*  and  thus  x*  G  A"*  are  optimal 
solutions  for  the  TN([f/°,  0]T,  [w,  1]T).  □ 

This  concludes  the  development  of  proposed  exact  methods  for  hireling  ^-efficient  deci¬ 
sions.  In  order  to  examine  their  practical  applicability,  the  following  section  discusses  several 
results  obtained  from  implementing  algorithms  for  the  constrained  objective,  the  reference 
point  and  the  set  lifting  approach. 

4  Numerical  examples 

We  demonstrate  the  above  methods  for  the  engineering  problem  of  designing  a  four  bar 
plane  truss  as  given  in  Figure  4,  taken  from  the  treatment  in  Koski  (1988),  and  compare  the 
results  to  those  obtained  from  applying  the  approximate  procedure  previously  proposed  in 
Engau  and  Wiecek  (2005b).  While  the  problem  implementation  for  that  paper  adopted  the 
formulation  given  in  Cocllo  Cocllo  (2001)  who  refers  back  to  Stabler  and  Dauer  (1992)  but 
introduces  a  slight  (and  probably  unintentional)  deviation  from  the  original  formulation,  we 
corrected  the  previous  implementation  for  all  following  results  accordingly. 
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The  problem  is  formulated  as  a  biobjective  program  with  the  two  conflicting  objectives 
of  minimizing  both  the  volume  V  of  the  truss  (/i)  and  the  displacement  A  of  the  joint  (f2), 
subject  to  given  physical  restrictions  on  the  feasible  cross  sectional  areas  x±,  x2,  x3  and  x4 
of  the  four  bars.  The  stress  on  the  truss  is  caused  by  three  loading  forces  of  magnitudes  F 
and  2 F  as  depicted  in  Figure  4.  The  length  L  of  each  bar,  the  Young’s  modulus  of  elasticity 
E  and  the  only  nonzero  stress  component  a  are  modeled  as  constants.  The  mathematical 
formulation  of  this  problem  is  given  as 


minimize 


fi(x)  =  L(  2xi  +  V2x2  +  \[2x3  +  x4), 
FL  (  2  2y/2  2^2  2 


h{x)  =  —  —  + 


E 


X\ 


X2 


X3 


H - 

£4 


subject  to  ( F jo )  <  x4  <  3 (F/a), 

V2(F/a)  <x2<  3 (F/a), 

\/2 (F/a)  <  x3  <  3(F / a), 


(F/a)  <  x4  <  3 (F/a). 


with  constant  parameters  F  —  10  kN,  E  —  2x  105  kN/cm2,  L  =  200  cm  and  a  =  10  kN/cm2. 

The  sets  of  attainable  outcomes  for  this  problem  and  e-Pareto  outcomes  as  identified 
by  the  approximate  weighted-sum  method  are  now  depicted  in  Figure  5,  where  a  detailed 
description  regarding  the  computation  of  the  latter  is  provided  in  the  original  paper. 

Observe  that  a  reasonable  relaxation  from  Pareto  to  e-Pareto  outcomes  does  only  occur  in 
the  middle  region  of  the  Pareto  curve  and  that,  in  general,  the  number  of  distinct  generated 
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Figure  5:  Sets  of  attainable,  Pareto  and  e-Pareto  outcomes  (approximate  method) 


outcomes  is  much  greater  for  the  Pareto  than  for  the  e-Pareto  outcomes.  Moreover,  many 
Pareto  outcomes  are  also  identified  when  solving  the  problem  with  specified  tolerances,  and 
it  turns  out  that  a  further  increase  of  the  specified  tolerances  for  truss  volume  and  joint 
displacement  yields  an  even  further  reduction  in  the  number  of  identified  e-Pareto  outcomes. 
As  mentioned  before,  we  believe  that  parts  of  these  methodological  issues  might  result  from 
the  difficulty  of  deciding  on  e-optimality  when  solving  the  SOP  due  to  insufficient  knowledge 
of  accurate  lower  bounds  on  the  optimal  solution.  Consequently,  although  the  obtained 
set  does  include  several  points  that  relax  efficiency  in  the  desired  manner,  in  general  the 
proposed  approximate  methods  do  not  provide  a  satisfying  representation  of  the  set  of  e- 
Pareto  outcomes  and  thus  reinforce  the  necessity  for  exact  methods  as  developed  in  this 
paper. 

Figures  6  and  7  show  the  results  obtained  when  applying  the  constrained  objective  and 
the  reference  point  method,  respectively,  for  the  four  different  tolerances  e1  =  (200,  0.002)r, 
e2  =  (500,  0.005)t,  e3  =  (800,  0.008)t  and  e4  =  (1000,  0.01)T.  As  intended,  both  approaches 
yield  a  reasonable  set  of  fully  relaxed  e-outcomes  for  each  of  these  four  different  values. 

Recall  that  the  formulation  of  the  constrained  objective  approach  guarantees  that  the 
identified  solutions  in  decision  or  outcome  set  are  truly  e-efficient  or  e-Pareto,  respectively, 
while  targeting  at  a  fully  relaxed  e-efficient  solution.  The  reference  point  method,  on  the 
other  hand,  merely  attempts  to  get  as  close  as  possible  to  some  e-Pareto  outcome,  thereby 
possibly  exceeding  those  e-constraints  imposed  by  the  constrained  objective  approach.  This 
observation  explains  the  different  behavior  of  these  two  methods  in  the  lower  right  regions 
of  Figures  6  and  7. 
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Figure  6:  Results  using  the  constrained  objective  approach  (exact  method) 


Another  difference  between  the  constrained  objective  and  the  reference  point  approach  is 
that  the  first  remains  a  multiple  objective  program  and  thus  requires  a  previous  scalarization 
before  subsequent  optimization.  Based  upon  the  scalarization  function  used,  this  might  lead 
to  additional  numerical  difficulties  possibly  explaining  the  few  outliers  not  occurring  when 
applying  the  reference  point  method.  In  general,  however,  both  approaches  resolve  the  exist¬ 
ing  drawbacks  of  the  approximate  methods  and  yield  an  almost  complete  representation  of 
the  set  of  e-outcomes  or  fully  relaxed  e-Pareto  outcomes  for  various  choices  of  the  tolerances 
e  over  the  whole  attainable  set. 

The  results  shown  in  Figure  8  are  obtained  from  implementing  the  weighted-sum  set 
lifting  approach  and  reveals  that  this  method,  although  conceptually  appealing,  is  subject 
to  various  numerical  issues  and  difficulties.  In  particular,  by  potentially  allowing  a  division 
by  zero,  this  formulation  becomes  very  sensitive  and  hence  does  not  seem  suitable  for  actual 
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Figure  7:  Results  using  the  reference  point  approach  (exact  method) 


implementation. 


5  Conclusions 

In  many  real  world  situations,  multiobjective  programs  can  only  be  solved  in  approxima¬ 
tion.  Since  the  set  of  efficient  decisions  commonly  consists  of  infinitely  many  points,  only 
an  approximate  representation  can  be  found  in  form  of  a  finite  set  of  points.  Furthermore, 
these  points  themselves  may  only  be  approximate  to  efficient  decisions  due  to  computational 
restrictions,  and  thus  the  importance  of  good  approximation  methods  in  multiobjective  pro¬ 
gramming  is  well  recognized,  see  Ruzika  and  Wiecek  (2005).  Nevertheless,  searching  the 
literature  suggests  that  applicable  methods  to  actually  identify  approximate  solutions  have 
not  been  well  developed. 
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Figure  8:  Results  using  the  set  lifting  approach  (exact  method) 


Following  preliminary  work  presented  in  Engau  and  Wiecek  (2005b)  which  also  empha¬ 
sizes  the  relevance  of  suboptimal  solutions  for  the  assessment  of  sensitivity  and  practical 
decision  making  situations  dealing  with  collections  of  more  than  just  one  multiple  objective 
program,  this  paper  continues  the  development  of  generating  methods  for  e-efficient  solutions 
and  gives  theoretical  proofs  for  their  capability  of  hireling  particularly  those  solutions  achiev¬ 
ing  the  full  relaxation  specified  by  the  tolerance  vector  e.  Application  of  these  methods  to 
an  engineering  design  problem  shows  that  various  drawbacks  of  those  method  derived  earlier 
are  resolved  and  therefore  leads  to  the  conclusion  that  the  new  methods  are  also  practically 
suited  for  that  purpose. 

We  believe  that  the  formulation  of  other  methods  for  generating  e-efficient  solutions  is 
only  limited  by  personal  creativity  and  thus  encourage  interested  researchers  to  propose  new 
approaches  towards  this  goal.  Our  personal  future  interest  lies  in  further  application  of  the 
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concept  of  e-efficiency  for  different  decision  making  situations  and  new  solution  methodolo¬ 
gies  for  multiple  objective  programming  problems. 
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